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Abstract 



An analysis of Deeply Virtual Compton Scattering (DVCS) is made within 
the colour dipole model. We compare and contrast two models for the dipole 
cross-section which have been successful in describing structure function data. 
Both models agree with the available cross section data on DVCS from HERA. 
We give predictions for various azimuthal angle asymmetries in HERA kine- 
matics and for the DVCS cross section in the THERA region. 



1 Introduction 



In this paper we explore the predictions of the colour dipole model for high energy 
deeply virtual Compton scattering (DVCS): 

l*(q)+p(P)^l(q')+p(P'), (1) 

where the first photon has spacelike virtuality q 2 = —Q 2 > 0, but the second photon 
is real (q' 2 = 0), and hence transversely polarized. DVCS is a particular example of 
a diffract ive process 

J*(q)+p(P)->X+p(P'), (2) 

in which the diffractively-produced system, X, is separated by a rapidity gap from the 
elastically-scattered proton (at least for high photon-proton centre-of-mass energies, 
W, i.e. W 2 = (q + P) 2 > Q 2 ,M],). The first HERA data on this process is now 
available |, |]. 

The colour dipole model of diffraction |j| provides a simple unified picture of such 
diffractive processes and enables "hard" and "soft" physics to be incorporated in a 
single dynamical framework. At high energies, in the proton's rest frame, the virtual 
photon fluctuates into a hadronic system (the simplest of which is a qq dipole) a 
long distance upstream of the target proton. The formation time of this hadronic 
system, and of the subsequent formation of the hadronic final state, is much longer 
than the interaction time with the target. It is this observation that leads to the 
main (plausible) assumption of the colour dipole model, i.e. that the interaction of a 
given fluctuation with the target is independent of how it is formed, and is therefore 
universal. It leads to the following generic factorization of the amplitudes of high 
energy diffractive processes: 

A(>y + p-+Xp) = J ^ n a^° ut (3) 

where a is the interaction cross section of a given configuration with the target and 
the integral runs over the phase space describing the incoming and outgoing hadronic 
systems. For the case of dipole scattering, one must integrate over dipole configura- 
tions (longitudinal momentum fractions and transverse sizes). We know of no formal 
proof of this type of high-energy factorization, whether applied to dipoles or more 
complicated configurations. Nevertheless, within this common framework there are 
many different formulations for the interaction cross section a [[3]1-|[11|, which have 
been applied with varying degrees of success^]. Here we consider two particular dipole 
models || ||| which have both been successful in describing structure function data, 
but which at first sight differ quite drastically in their structure and implications; and 
compare their predictions for DVCS. 

*For a recent overview, see fl2l. 



An additional assumption of most dipole models of diffraction is that the scattering 
with the target is diagonal with respect to the appropriate variables (i.e. transverse 
sizes, momentum fractions and polarizations are unchanged by the interaction). For 
the case of DVCS this implies that the incoming photon must be transversely polarized 
in order to respect s-channel helicity conservation. 

DVCS is a good probe of the transition between soft and hard regimes in the dipole 
model for two reasons. Firstly, the transverse photon wave function can select large 
dipoles, even for large Q 2 , and certainly for the Q 2 range 2 < Q 2 < 20 GeV 2 for which 
data is now available [fj]. Secondly, because the final photon is real, DVCS is more 
sensitive to large dipoles than DIS at the same Q 2 , as we shall illustrate quantitatively 
in Section 3. In addition, for Q 2 — > 0, the process reduces to real Compton scattering 
and the cross-section can be reliably inferred from real photo-absorption data, where 
soft physics dominates. 

We stress the potential importance of well-founded dipole descriptions in providing 
reliable starting points for exploiting DGLAP evolution properties at "large" Q 2 . 
From the theoretical point of view, DVCS is the best understood of all exclusive 
diffr active processes, essentially because the X system is just a real photon. Indeed, 
a perturbative QCD factorisation theorem has been explicitly proven as Q 2 — > oo 
|13| which enables the QCD amplitude to be described by a convolution in momen- 



tum fraction of generalised (or skewed) parton distributions [|TJJ (GPDs) with hard 
coefficient functions. GPDs correspond to Fourier transforms of operator products 
evaluated between proton states of unequal momenta (cf. eq.(|])). They are there- 
fore generalizations of the familiar parton distributions of deep inelastic scattering, 
and like them satisfy perturbative evolution equations [15, [16|, |T7|, 18] which enable 



them to be evaluated at all Q 2 in terms of an assumed input at some appropriate 
Q 2 = Qq. In practice, to compare with experimental results at finite Q 2 one must 
establish a regime in Q 2 in which the higher twist corrections (see e.g. |ID| , p]|] ) to 
this leading twist result are numerically unimportant. This is a very difficult task in 
general but early indications are that the minimum Q 2 values defining this regime 
are considerably higher than in inclusive cross sections, for which values as low as 
Qq = 1 GeV 2 , or even lower, have been used. Since the contributions from different 
transverse sizes are manifest in the dipole model one may realistically hope to gain 
insight into this question by investigating DVCS in the dipole framework. As such 
a dipole analysis of DVCS provides a complementary description to the formal QCD 
analysis, applicable at "large" Q 2 . Any insight gained regarding the mixture of soft 
and hard physics within the dipole model framework, can also be employed in those 
processes for which factorization theorems have not been proven. 



Frankfurt, Freund and Strikman |2T[ have given a leading order QCD analysis of 
DVCS. The resulting predictions for the DVCS amplitude at t — are in agreement 
with the recent HI measurements M of the total DVCS cross-section, assuming an 



exponential t- dependence with a reasonable value of the slope parameter. The GPDs 
are evolved from an input value Qq = 2.6 GeV 2 , where the input GPDs are obtained 
by estimating their ratio to "ordinary" parton distribution functions (PDFs) using a 
simple aligned jet modef]. While this provides a reasonable first estimate, it is clearly 
subject to uncertainties which will become important when more accurate data are 
available. 

Recently NLO QCD analyses of DVCS have been completed fli~8 , |23| , p4] , g5j which use 



as input GPDs Radyushkin's model [26] based on Double Distributions proportional 
to PDFs (which automatically impose the correct symmetry properties in the so- 
called "ERBL region"). The colour dipole model offers a means of estimating these 
distributions at the input scale in the DGLAP region, in a complementary and well- 
founded framework, which can accurately describe both virtual Compton scattering 
and other closely related data over a wide range of Q 2 . This is possible because at 
leading- log accuracy in Q 2 , the amplitude is approximately equal to the GPD, at a 
particular point. In this paper we compute predictions for the cross section of eq. ([[]), 
and for various azimuthal angle asymmetries for the associated lepton process |p7| 
which are sensitive to both the real and imaginary parts of the DVCS amplitude. 

The structure of the paper is as follows: in section ^| we summarize and compare 
our two dipole models in the context of deep inelastic scattering; we then discuss 
their application to virtual Compton scattering process in section ||]; compute various 
observables in section |] and summarize our results and conclusions in section |5|. 



2 The colour dipole model 

Singly dissociative diffractive jp processes (cf. eq.(|])) are conveniently described in 
the rest frame of the hadron, in which the incoming photon dissociates into a qq pair 
a long distance, typically of order of the "coherence length" 1/Mx, from the target 
proton. Assuming that the resulting partonic/hadronic state evolves slowly compared 
to the timescale of interaction with the proton or nuclear target, it can be regarded 
as frozen during the interaction. In the colour dipole model, the dominant states are 
assumed to be qq states of given transverse size, dr- Specifically 

|7r) = fdz d 2 d T tf(z, d T , Q 2 ) \z, d T , ) + . . . , (4) 

where z is the fraction of light cone energy carried by the quark and ip?(z, dx, Q 2 ) is 
the light cone wave function of the photon of polarization r = T,L. Assuming that 
these dominant states are scattering eigenstates (i.e. that z, dx and the quark helic- 
ities, which are left implicit in the above equation, remain unchanged in diffractive 



^For further discussion of this approximation, see p2[ and the original paper El 



scattering) the elastic scattering amplitude for j*P ~^ 1*V is specified by Fig. |T]. This 
leads via the optical theorem to 



a 



i p 



Jdzd 2 d T \^l L (z,d T} Q 2 )\ 2 a(s* } d T} z) , (5) 



for the 7*p total cross-section in deep inelastic scattering, where a(s*, d?, z) is the 
total cross-section for scattering dipoles of specified (z, d T ) which do not change in 
the interaction (the second line then follows from orthogonality and the variable s* 
will be specified shortly). 



Figure 1: The colour dipole model for the elastic process 7*p — > 7*p(DIS), and virtual 
Compton scattering 7*p — > jp. 

The dipole cross-section is usually assumed to be flavour independent and "geometric" , 
i.e. independent of z. Beyond this the models fall into two main classes. 

In the first, the dipole cross-section is assumed to depend solely on the properties of 
the dipole-proton system itself, implying the choice s* = W 2 . Other singly diffractive 
photo-processes involve exactly the same dipole cross-section, but different wavefunc- 
tion factors depending on the final state, as we shall see below for virtual Compton 
scattering. 

The second type is more closely connected with hard perturbative QCD predictions 
for the interaction cross section, which for small x and high Q 2 involve two gluons 
being exchanged. For hard scattering, small dipoles are connected via two parton 
lines to the proton. The interaction cross section then depends on the momentum 
fractions of the proton carried by the parton lines, i.e. s* = xbj = Q 2 /W 2 , or 
s* = x' ~ xbj- In this case the dipole cross-section must be slightly modified when 
applying it to different processes, since GPDs of the appropriate kinematics must be 
used, as discussed below. 

In the rest of this section, we shall briefly summarize the properties of one model of 
each type, after first considering the other main ingredient in eq.([5[), i.e. the photon 
wavef unction. 



2.1 The photon wavefunction 

Because the proton structure function, F 2 , is predominantly transverse, both small 
and large dipoles contribute significantly to F 2 over a wide range of Q 2 , where "large" 
means transverse sizes of order d? « 1 fm. With this caveat, it is none the less useful 
to consider two qualitatively different regimes. 

For small dipoles, it is reasonable to assume "QED wavef unctions," ^ TL = 4>tl 
calculated^ from the usual QED vertex —icy 11 . Explicitly 

nr 
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\^ L (z,d T ,Q 2 )\ 2 = -a e . m .Ve 2 QV(l-z) 2 K V T ) (6) 

q=l 

rif 

\^ T (z,d T ,Q 2 )\ 2 = l- 2 a e . m .J2e 2 q {[z 2 + (l^z) 2 }e 2 K 2 (ed T )+m 2 q K 2 (ed T )}(7) 



q=l 



where 



e 2 = z(l-z)Q 2 + m 



<i ■ 



K and K\ are modified Bessel functions and the sum is over quark flavours. Fur- 
thermore, for the large Q 2 values where small dipoles dominate, these wavefunctions 
become insensitive to the quark mass. In this regime, the wavefunctions are essentially 
known. 

For small Q 2 — 4m 2 the QED wavefunctions become sensitive to the squared quark 



q 



mass m . At the same time, large dipoles, for which one would expect significant 



confinement corrections, become very important. The wavefunction is clearly model 
dependent in this region. 

Both the models discussed in this paper assume that the dipole cross-section becomes 
"hadron-like" for dr~ 1 fm, with an energy dependence characteristic of the "soft 
Pomeron" . In choosing the wavefunction in this region, FKS || were motivated by the 



work of Frankfurt, Guzey and Strikman [29|. These authors analyzed the distribution 



of scattering eigenstates in a non-diagonal generalized vector dominance model [30 



which provides a good description of the soft Pomeron contribution to the nucleon 



structure function F 2 on both protons and nuclei |31|. They found a distribution 
of states which was qualitatively similar to that obtained in a colour dipole model 
with a perturbative wavefunction, but with an enhanced contribution from dipole 
cross-sections of hadronic size. In the light of this, FKS chose m 2 = 0.08 GeV 2 
corresponding roughly to a constituent mass for the light quark case; and modified 
the QED wavefunction by multiplying by an adjustable Gaussian enhancement factor: 

\^T,L(^dTiQ 2 )\ 2 = \^° T Az,d T ,Q 2 )\ 2 f(d T ) (8) 



■i-For an explicit derivation, see Appendix A of 
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Figure 2: The weight function H(cIt) for different Q 2 corresponding to the photon 
wavefunction obtained by FKS || with m 2 = 0.08 GeV 2 . The peak at low Q 2 
represents the modification to the perturbative photon wave function at large cIt- 



0.5 - 




d T [fm] 



Figure 3: The weight function H{dr) for various Q 2 corresponding to photon wave- 
functions of perturbative form for m 2 = 0.08 GeV 2 (solid lines) and m 2 = 0.02 GeV 2 
(dotted lines). 



where 

f{dT) ~ l+BeM-c 2 R 2 ) ' ( } 

This form enables the width and height of the enhancement to be controlled inde- 
pendently while keeping a factor of close to unity at small dr- The effect of this is 
conveniently summarized by integrating out the angular and z dependence in eq.@ 
to give 

oli = Jdzd 2 d T (\Mz,d T )\ 2 + \Mz,dT)\ 2 )t(s,d T ) 

12 f 

= — oi em , \ dd T H(d T )a(s,d T ) . (10) 

The resulting behaviour of H{dr) for the final parameter values (see below) is shown 
in Fig. 0: as can be seen, the enhancement is important for very low Q 2 , but de- 
creases rapidly as Q 2 increases. Other authors do not in general include an explicit 
enhancement factor, but achieve a similar effect, at least for Q 2 > 1 GeV 2 , by varying 
the quark mass. Choosing a smaller quark mass increases the wavefunction at all 
large dx, as illustrated in Fig. |3|. Golec-Biernat and Wusthoff ||, for example, used 
m 2 = 0.02 GeV 2 , comparable with the pion mass squared. 

For our second model, the MFGS model ||, this question is less important, since 
results are only presented for Q 2 > 1 GeV 2 ; and results are presented for m 2 = 
0.08 GeV 2 without a confinement factor of eq.(|^). In both models, a charmed quark 
contribution has also been included, which only differs from the up quark contribution 
by the quark mass m 2 = 1.4 GeV 2 . 

Given this uncertainty in the wavefunction, it is clear from eq.(0) that the dipole 
cross-section at large dr ~ 1 fm cannot be inferred, even in principle, from structure 
function and real photo-absorption data alone. Other information must also be used. 



2.2 The FKS model 

This model belongs to the class in which the dipole cross-section is assumed to 
depend solely on the properties of the dipole-proton system itself, implying the choice 
s* = W 2 independent of the virtuality of the incoming (or outgoing) photon. The idea 
was then to extract information on the dipole cross-section by assuming a reasonable 
but flexible parametric form to fit structure function and real photoabsorption data 
in the diffractive region x < 0.01 for < Q 2 < 60 GeV 2 . This was implemented by 
assuming a sum of two terms 

a(W 2 , d T ) = a soh (W 2 , d T ) + a haid (W 2 , d T ) , (11) 



each with a Regge type energy dependence on the dimensionless energy variable 
d 2 T W 2 : 

° S oit(W\d T ) = a s (l - IT ^)(4^ 2 ) As (12) 

^hard^ dr) = (a$ ' d% + a%d%) exp(-u H d T )(d 2 T W 2 ) XH (13) 

These functions were chosen^ so that for small dipoles the hard term dominates yield- 
ing a behaviour 

a -»■ a%d 2 T (d 2 T W 2 ) XH d T ^0 

in accordance with colour transparency ideas. For large dipoles the soft term domi- 
nates with a hadron-like behaviour 

a ss a s (d 2 T W 2 ) Xs d T « 1 fm. 

The values As ~ 0.06, Xh ~ 0.44 resulting from the fit are characteristic of the 
soft and hard Pomeron respectively, but the fits could be obtained for a range of 
values for the parameter a^ because of the uncertainty in the photon wavefunction at 
large dx discussed above. This ambiguity can be resolved by using the same dipole 
cross-section to calculate the structure function, F 2 \x, Q 2 , M|), for diffractive deep 
inelastic scattering (DDIS) 

7*+p->X + p (14) 

and a subsequent paper M showed that good agreement was found for af ~ 30 GeV~ 2 . 
The parameter values for this fit are given in Table 1 and the resulting behaviour of 
the dipole cross-section as a function of d? is shown in Fig. f| for three energies, 
including W = 75 GeV corresponding to the mean energy of the virtual Compton 
scattering data |2j to be discussed below. 

One feature of the FKS model is that, in its present form, it does not include "gluon 
saturation" or "unitarity corrections" which are expected to eventually damp the 
rapid rise with energy of the dipole cross-section for small dipoles. Its success implies 
that such effects are not necessarily required in the HERA region. However, by 
examining the predicted behaviour of the dipole cross-section, the authors have argued 
H that saturation effects will begin to play a role at the top of the HERA range and 
will rapidly become important above it. This is illustrated in Fig. [|, which shows that 
the rapidly rising dipole cross-section at small dr, where the hard term dominates, 
reaches hadronic sizes at the top of the HERA range W ~ 300 GeV. 

Finally, for completeness, in Fig. |5] we compare the FKS dipole cross-section with 



that used in a recent analysis of DVCS by Donnachie and Dosch [O]. This model 

§In B a more complicated parametric form was used, but this simpler parametric form gives 
very similar results. 
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Table 1 



As 


0.06 ±0.01 


Xh 


0.44 ±0.01 


c*o 


30.0 (fixed) 


of 


0.072 ±0.010 


af 


0.027 ±0.007 


of 


1.89 ±0.03 






vh 


3.27 ±0.01 


B 


7.05 ±0.08 


c 2 


0.20 (fixed) 


R 


6.84 ±0.02 






ml 


0.08 (fixed) 


2 


1.4 (fixed) 



Table 1: Parameters for the FKS model specified by eqs.(||, |Tl|), in appropriate GeV 
based units throughout. 
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Figure 4: The FKS dipole cross section at W = 10, 75, 300 GeV (dashed, solid and 
dotted lines, respectively). 
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Figure 5: The Donnachie-Dosch dipole cross section ||TT| at W = 75 GeV (dashed 
line) and W = 300 GeV (dotted line). The FKS dipole cross-section at W = 75 GeV 
(solid line) is included for comparison. 



associates dipoles whose sizes are less than (or greater than) a certain critical size 
with a fixed power energy dependence corresponding to the hard (or soft) Pomeron 
respectively^} For dipoles whose transverse size is less than about 0.8 fm, the resulting 
behaviour is not unlike the FKS dipole cross-section at W = 75 GeV, which is the 
mean energy of the data, and the model gives a good account of HI data in this 
region. However at higher energies the model develops a rather artificial dependence 
on the dipole size: as can be seen from the dotted curve in Fig. |5| for W = 300 GeV, 
the dipole cross section not only does not increase monotonically with dipole size, 
but develops a large discontinuity at the matching point (d± ~ 0.3 fm). The other 
obvious difference is that the cross-section is much larger than the FKS cross-section 
at large dr- However for DIS this can be compensated by differences in the photon 
wave-function, which is uncertain in this region, as noted at the end of section 2.1. 
FKS resolved this ambiguity by considering the DDIS reaction of eq. flT4|) , which 
is more sensitive to large dipoles and involves a strikingly different combination of 
wavefunction and cross-section. It would be interesting to see the predictions of the 
Donnachie-Dosch model for this reaction. 



2.3 The MFGS model 



This model |], [TIJ is directly based on the known behaviour of hard small- x QCD 
processes, i.e. that they are driven by the gluon distribution at small x. Using the 
phenomena of colour transparency, it directly relates the dipole cross-section at small 
dr to leading order (LO) gluon distributions at large Q 2 . To leading order in InQ 2 , 

' In contrast the FKS model has both hard and soft components for all sizes, but in smoothly 
varying amounts. 
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and within the small x limit, the total photon-proton cross-sections are given by 
expressions of the form eq.(f|) with a QED wavefunction and a dipole cross-sectionf] 



o 



ir 2 d 2 

T (r\2\ i i r\2\ 



pQCD 



[x, d T ) = —^a s (Q 2 )xg(x', Q 2 ) , (15) 



where xg(x',Q 2 ) is the LO gluon distribution of the proton. At leading log it is 
sufficient to choose x' = x and Q 2 = Q 2 . However in the MFGS model (specified fully 
in ||) an attempt to go beyond leading log was made by introducing <i T -dependence 
into the scales, x' and Q 2 . For the four momentum scale Q 2 , the phenomeno logical 
relation 



«• = ;? (16) 

was assumed where A =<d%> Q 2 . A theoretical procedure for defining <d%> using 
the integral in dx for Fl in |33j gave a value of A ~ 10. This value was used in ||, 
but it was later discovered that the inclusive cross sections are rather insensitive to 
its precise value in the range A = 4 — 15 and that the lower value of A = 4 appears 



to be favoured by the J/^-photoproduction data [fLOX We adopt this lower value in 
what follows. The momentum fraction required to create a quark-antiquark pair of 
mass M 2 q = (k 2 ^ + m 2 )/(z(l — z)) is 

M 2 +0 2 
X Q 2 + W 2 ' {U) 

since kr is Fourier conjugate to dx the following relationship was adopted for a 
dipole of given transverse size dr'- 



x = x 



1 + 0.75- 



d\{Q 2 + Am 2 ) 
The dipole cross-section was then evaluated using the CTEQ4L gluon distributions 



for a region d? < d Tc in which eg. (jig) is appropriate, where the boundary d Tj( 



is specified below. For large dipoles dx > <^t> = 0.65 fm, the form 

*Wr > 4..,) = *(*P)^r {f)™ • (19) 

with an x-dependence characteristic of soft Pomeron exchange, is used, where Xq = 
0.01 and the value at dr = dr,n , x = xo is matched to the pion-proton total cross- 
section, a(nP) = 24 mb. In the intermediate region g?t, c < dr < c^t.tt, the dipole 
cross-section was linearly interpolated between the boundary values at g?t, c and g^tt- 

For moderate x, the point d?^ is set by the boundary of the perturbative region in 
dr- dT, c = dxfl — vVQo = 0.246 fm, for Q = 1.6 GeV which is the starting scale 
"This formula is implicit in most perturbative two gluon models. For an explicit derivation see 
©■ 
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Figure 6: The MFGS dipole at W = 10, 75, 300 GeV (dotted, solid and dashed lines, 
respectively) at fixed Q 2 = 1 GeV 2 , corresponding to x ~ 10~ 2 ,2.10~ 4 and 10 -5 , 
respectively. 



for the CTEQ4L partons (below this scale they are not defined). The resulting dipole 
cross-section for Q 2 = 1 GeV 2 is shown in Fig. || for various values of W. 

At small enough x, as a result of the rising gluon density, the small dipole cross 
section increases faster than the large dipole one and soon reaches hadronic size (tens 
of mb). This threatens to spoil the monotonic increase of a with d?- To prevent this 
the MFGS model implements taming corrections that guarantee that the small dipole 
cross section cannot reach more than half its value at rfr j7r . This constraint implies a 
dx^c that shifts to increasingly small d?- This correction is not crucial in the HERA 
region for A = 4, but does becomes important above it. 

The parameters of the model are not adjusted to fit data, but nonetheless good semi- 
quantitative accounts of the deep inelastic scattering |§ and J/ip photoproduction 
data [10] were obtained. For exclusive diffractive processes, such as vector meson 



production or DVCS, it is necessary to include GPDs, parameterized in terms of 
skewedness, 5, x' and Q 2 , rather than the ordinary ones used in eq.([H|). For DVCS, 
5 = xbj = Q 2 /W 2 . To implement leading order GPDs, we adapted the skewed 
evolution package developed by Freund and Guzey PS] , using the CTEQ4L gluon 
distributions [0 as input to the LO skewed evolution. 



The model focuses on small dipoles, and while the behaviour for large dipoles is a sen- 
sible guess, no detailed attempt to resolve the intrinsic ambiguity in the wavefunction 
and cross-section for large dipoles discussed earlier has been made. In what follows 
we shall restrict this model to Q 2 > 1 GeV 2 , where one is relatively insensitive to this 
region due to the smallness of the wavefunction. 
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3 Virtual Compton scattering 



In the colour dipole model, virtual Compton scattering is again given by Fig. [I], but 
with a real photon in the final state, leading to 

Im A DYCS (W 2 , Q 2 , t = 0) = J dz d 2 d T ip^(z, d T , Q 2 ) a(s*, z, d T ) ip T (z, d T , 0) , (20) 

for the imaginary part of the DVCS amplitude at zero momentum transfer. Thus, 
our dipole models yield no-free-parameter predictions for this process. In this section 
we compare the predicted behaviours of the amplitude in the two models, leaving the 
comparison between predictions and experiment to section [5]. 




d T [fm] 



Figure 7: Profile in transverse dipole size for Q 2 = 1 GeV 2 and W = 75 GeV, 
employing the FKS (solid line) and MFGS (dashed line) models for the dipole cross 
section. 
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Figure 8: Profile in transverse dipole size for Q 2 = 10 GeV 2 and W = 75 GeV, 
employing the FKS (solid line) and MFGS (dashed line) models for the dipole cross 
section. 
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We start by comparing the contributions to the amplitude arising from dipoles of 
different size. To do this we perform the angular and z integrations to rewrite eq.fpOp 
in the form 

Im A{W 2 , Q 2 , t = 0) = 2vr f dd T p{d T , s, Q 2 ) , (21) 

where the profile function 

p(d T , s*, Q 2 ) = J dzd T tP* T (z, d T , Q 2 ) &{s*, z, d T ) ip T {z, d T , 0) , (22) 

gives the relative contributions arising from dipoles of different size dr- The results 
are shown for both models at the mean energy, W = 75 GeV, of the HI data in 
Figs. [7], |8|. As Q 2 increases, the profile shifts to smaller dx- The FKS model has a 
larger contribution from large dipoles than the MFGS model, although the forward 
amplitudes, obtained by integrating over all transverse sizes, are similar over a wide 
range of W and Q 2 , as we shall see. 

It is also interesting to compare eq.(p2|) to the corresponding profile function 

p T (d T , s*, Q 2 ) = dz d T iPt(z, dr, Q 2 ) a(s*, z, d T ) ipr(z, d T , Q 2 ) , (23) 



for the forward Compton scattering amplitude in which both transverse photons have 
the same Q 2 , which is related by the optical theorem to the transverse cross-section 
in DIS. The characteristic behaviour differences observed between DIS and DVCS are 
illustrated in Figs. [5] and W, using the FKS dipole model. For Q 2 = the two profiles 



are obviously identical, and for all Q 2 they become identical for small dr, since to 
leading order in 1/dr, the transverse wavefunction 



mz,d T ,Q 2 )\ 2 -> ^a e . m .X> f 2 + % Z? (24) 

9=1 T 

independent of Q 2 . However, as can be seen, when Q 2 increases the large dipoles 
are less suppressed in the DVCS case than in the DIS case, so that the former is the 
"softer" process at any given Q 2 . 

Returning to DVCS, the imaginary part of the amplitude is trivially obtained by 
integrating eq. (|2l|) over all dr- In the FKS model we have a sum of two Regge 
contributions, and the real part is easily computed from the corresponding signature 
factors; for the MFGS model, results were obtained using dispersion relations, as in 
||10||. The results for real and imaginary parts, using both models, are plotted as a 



function of W in Fig. 11 at the mean Q of the HI data, and the ratio is plotted in 



Fig. [12]. We extend the energy range out to the THERA range (see e.g. P5| ) and 
one can clearly see that the FKS model has a steeper energy dependence at very high 
energies. 
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Figure 9: Comparison of the profile functions of eqs.(|22|, |23|) for DVCS (solid line) 
and transverse DIS (dot-dash line) respectively at Q 2 = 1 GeV 2 and W = 75 GeV, 
employing the FKS model for the dipole cross section. 
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Figure 10: Comparison of the profile functions of eqs.(|22|, p3|) for DVCS (solid line) 
and for transverse DIS (dot-dash line) at Q 2 = 10 GeV 2 and W = 75 GeV, employing 
the FKS model for the dipole cross section. 
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Figure 11: Real and imaginary parts of the DVCS amplitude for the FKS (solid lines) 
and MFGS (dashed lines) dipole models for Q 2 = 4.5 GeV 2 . 
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Figure 12: The ratio, /3, of the real to imaginary parts of the DVCS forward amplitude 
at Q 2 = 4.5 GeV 2 : FKS (solid line); MFGS (dashed line); and FKS at Q 2 = GeV 2 
(dotted line). 
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Finally, in our introduction we noted that in |Hj], the ratio 



R = 



Im A{i*N -> 7*iV) t=0 
Ira A(j*N -* jN) t=0 



(25) 



of the imaginary parts of the forward amplitudes for DIS and DVCS was estimated at 
the input Q 2 = 2.6 GeV 2 using a simple aligned jet model, in order to infer the input 



generalized parton distributions for QCD evolution. Explicitly, this model gives [21 



R 



Q 2 



Q 2 + M 2 



\n-\l + Q 2 /M 2 ) 



(26) 



where Mq is estimated to be in the range 0.4 - 0.6 GeV 2 . The predictions of our 



models are compared with eq.(^) in Fig. [13], suggesting somewhat larger values at 
Q 2 = 2.6 GeV 2 , as can be seen. 
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Figure 13: The ratio of the forward amplitudes for DIS and DVCS (cf. eq.(p5|)) at 
W = 75 GeV: FKS (solid line); MFGS (dashed line); and the prediction of eq.([2T]) 
for Ml = 0.6 GeV 2 (dotted line) and M 2 = 0.4 GeV 2 (dot-dashed line), respectively. 



4 The lepton level process 



4.1 Definition of the DVCS cross section 



Virtual Compton scattering is accessed experimentally through the leptonic process: 

e ± (k) + p{p) -► e ± (k') + p{p') + 7 (g') (27) 

where the four momenta of the incoming and outgoing particle are given in brackets. 
As well as DVCS, the Bethe- Heftier process (BH), in which the photon is radiated by 
the initial or final state lepton, also contributes. On integrating over the azimuthal 



angle (denned below), the interference term between the two processes vanishes in 
the limit of large Q 2 , and the differential cross-section can be written as: 

d 2 a d 2 a DVCS d 2 a BH 



+ 



dydQ 2 dydQ 2 dydQ 2 

where y = (k - k') ■ P/{k ■ P) and Q 2 = -(k - k') 2 . For DVCS, Q 2 is the magnitude 
of the virtuality of the (spacelike) virtual photon and, in the proton's rest frame, y is 
the fraction of the incoming electron energy carried by the virtual photon. Neglecting 
the lepton and proton masses we have y ~ (Q 2 + W 2 )/S, where S — (k + p) 2 is the 
square of the lepton-proton centre-of-mass energy. The Bethe- Heftier contribution is 
essentially known in terms of the Dirac and Pauli form factors (see e.g. eqs.(18, 27) 
of | f27|| ) and can be easily calculated and subtracted from the total to leave 

d 2 DVCS 

w^^ |1+(1 - y)V(7 ' p ^ 7P) - (28) 

Making a trivial change of variable from y to W = \/(k — k' + P) 2 yields: 

(f DVCS 

W = ^ [1 + (1 " !/)V(7 *^^ ) ' (29) 

We now have a convenient form for comparing to the data on cr(7*P — > 'jp) which is 
binned in Q 2 and W. Assuming the usual exponential dependence in t = {p—p') 2 < 0, 
i.e. e Bt , the total 7* — p cross section is given by: 

°{l*P -»• IP) = 7^ -^ U=o , (30) 

where, with our definition for A DVCS , 

and j3 is the ratio of the real to the imaginary part of the forward virtual Compton 
scattering amplitude A at t — 0. 



4.2 Definitions of asymmetries 

A unique and attractive feature of DVCS is the interference with the Bethe-Heitler 
process which offers the rare chance to isolate both real and imaginary parts of the 
diffractive amplitude via azimuthal angle asymmetries [E7J. These asymmetries are 



conveniently discussed in a special frame P| with the proton at rest such that the 



direction of the vector q = k — k! defines the negative z-axis. Without loss of generality 
we can choose the incoming electron to have only a non-zero component along the 
positive x-axis in the transverse (x — y) plane. In this frame we have the following 
four-vectors: 

k= (k ,k sin 6 e ,0,k cos 9 e ), q = (q ,0,0,-\q 3 \) (32) 

P= (M, 0,0,0), P' = (P ', |P'|cos0sin0 H , |P'| sin0sin0 H , |P'| cos0 H ), (33) 
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where the angle of interest, 
hadron scattering planes. 



is the azimuthal angle between the lepton (x — z) and 



The motivation for using this frame is that the frame- dependent expression for the 
u— channel BH lepton propagator has a particularly simple Fourier expansion in the 
angle <fi. In [E7|| a slightly different frame is used and an explicit expression for the 
u— channel propagator is given, up to terms of order 0(1/Q 3 ) (cf. eq.(21) of |27j ): 



(k-q 



[i-y)Q 2 



(A + A 1 cos <p + A 2 cos 2<p + 



(34) 



where 



A 



A, 




(l-x)(l-2t min /t) 

(i-y) 



l-t 



mm 



/t)(l-x) 



y) 



(35) 



In the frame used here [19] the Fourier series terminates at cos0 (A n = 0, n > 2). 
We explicitly include factors of 1/(Aq + Aicos^) as appropriate in our numerical 
results. We use a code written for |23J which approximately implements eq.(18) and 
eqs. (24,27,30) of |27[ (the latter neglect terms of 0(1/Q), i.e. they use A = 1 and A 1 
= 0). The code includes the above expansion of the u— channel BH lepton propagator 
in our frame (taking its full and y-dependence into account, up to corrections of 
0(1/Q 3 )). The asymmetries of interest involve the quadruple differential cross section 
on the lepton level 

da Dvcs+BH = rfg (4) (ep~^ep7) 
dxdQ 2 dtdcf) 

In order to proceed it was necessary to convert our amplitude (cf. eq.(f20"D) to the 
dimensionless unpolarized helicity non-flip amplitude, Hi, appearing in eqs. (24,27,30) 
of [ ^7[j , since at small x and moderate t the contributions of the polarized, 7ii,£i, 
and unpolarized helicity-flip, S\, DVCS amplitudes are negligible. Assuming a simple 
exponential t— dependence on the amplitude level the conversion factor is: 



Hi = cxp( J Bt/2) 



W 2 



47T<X 



_ /idipok 



where W is in units of GeV, so the standard conversion factor fie 
is necessary to make Hi dimensionless. 



(36) 



0.389 GeV 2 mb 



Using the special frame defined above to specify the azimuthal angle, the asymmetries 



are defined as follows (see also p4|): 



The (unpolarized) azimuthal angle asymmetry (AAA), measured in the scat- 
tering of an unpolarized probe on an unpolarized target, is defined by 
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fj\ dcj>(da DVCS+BH - da BH ) - J 3 , 7 ; 72 dct>{da DVCS+BH - da BH ) 

AAA = — — (37) 

j^ d(j){da DVCS + BH - da BH ) 

where da BH is the pure BH term. The above approximation for the u— channel 
BH propagator leads to a non-trivial 0— dependence of the pure BH term. To 
directly access the DVCS amplitudes, via the interference term, we define AAA 
with this piece subtracted. Note that with this "subtracted definition" the 
magnitude of AAA may become greater than unity in certain regions. 

The single spin asymmetry (SSA), measured in the scattering of a longitudinally 
polarized probe on an unpolarized target, is defined by 



£ d<j)Aa DVCS+BH - Xf dcj>Aa DVCS+BH 
J^ d<j)(da DVCS+BH ^ + da DVCS + BH ^) 



SSA = ^^ J^T^ , (3 j 

pZty ii/i 7~it//^c i Dirt . i ni/pc i du i\ 



where Aa = do^ — da^ and j and J, signify that the lepton is polarized along or 
against its direction of motion, respectively. 

The charge asymmetry (CA) in the scattering of an unpolarized probe on an 
unpolarized target: 

f* /2 /2 d<j>Ad c a DVCS+BH - f 3 J^d^Ad c a DVCS+BH 

CA = —H ^ , (39) 

J^ d(j)(d+a DVCS + BH + d-a DVCS + BH ) 

where Ad c a = d + a — d~a corresponds to the difference of the scattering with 
a positron probe and an electron probe. 



In the small-x limit, required for the dipole approximation, and at large Q 2 , AAA 
and CA are directly proportional to the real part of the DVCS amplitude, and SSA 
to the imaginary part of the unpolarised amplitude (cf. eq. (30,40,43) of p7|)[**|. 



5 Results and conclusions 



The Q 2 and W dependencies of the DVCS total cross-section of eq. (ffDD , obtained 
from our two dipole models, are compared with the HI data in Figs. |T^, |H]. In 



these figures, the vertical errors include statistical and systematic errors added in 
quadrature, and the horizontal errors indicate the bin width of the data. In making 



**At large x one also gets a contribution from the polarised and helicity-fiip DVCS amplitudes, 
li\,£\,£i, which are negligible in the small x region (cf. eqs.(24, 30) of [p7|). 
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Figure 14: The energy dependence of the photon level DVCS cross section at fixed 
Q 2 = 4.5 GeV 2 : FKS (solid line); MFGS (dashed line). 
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Figure 15: The Q 2 -dependence of the photon level DVCS cross section at fixed W 
75 GeV: FKS (solid line); MFGS (dashed line). 
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Figure 16: Compendium of results on asymmetries for fixed x — 10 2 , at two values 
of Q = 2, 3 GeV, accessible in the HERA kinematic range. 
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Figure 17: Compendium of results on asymmetries for fixed x = 10~ 4 , at two values 
of Q = 2, 3 GeV, accessible in the HERA kinematic range. The results for SSA and 
AAA may be compared to Fig. 1 of 



and for CA to Fig. 4a of |7| 
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this comparison we have used the same value for the slope parameter, B = 7 GeV~ 2 , 
as that used by HI to take into account the resolution and acceptance of their detector. 
It should be emphasized that the uncertainty in B implies an associated uncertainty 
in the normalization of the predictions, which could well be weakly Q 2 dependent. In 
addition, Donnachie and Dosch |ll| have stressed that the real photon limit provides 



an important constraint, since the forward imaginary part can be inferred from fits to 
real photoabsorption data using the optical theorem, and their resulting value for the 
real photon cross section a^p — > jp) is also shown in Fig. |15p| . As can be seen, the 



consistency between the predictions and experiment is very encouraging for the basic 
features of the dipole models, despite the large statistical errors and the uncertainty 
in the slope parameter. Conversely, the close agreement between the predictions of 
the models out to energies well beyond the HERA region, despite the fact that one 
implements saturation effects approximately and one does not, means that the total 
cross-section is not a good discriminator between them. This is not unexpected, since 
DVCS is clearly a "softer" process than DIS, due to the presence of a real photon. 



In figs.(^q,PD we show results for the predicted asymmetries of eqs.(g7|, |38|, |39|) result- 
ing from both models to indicate the overall size and spread of predictions expected in 
the small- X region (x G [10~ 2 , 10~ 4 ]) of HERA kinematics. The AAA and SSA results 
are for a positron and we use HERA kinematics with proton energy E p = 920 GeV 
(i.e. S = AE e E p « 99, 000 GeV 2 ) to compute the value of y = (W 2 + Q 2 )/S. 

In conclusion, we observe that both models provide a good description of the avail- 
able DVCS cross section data, without further tuning. We give predictions for the 
cross section at higher energies. A measurement of the asymmetries would allow the 
predicted phase of DVCS amplitude for both models to be tested. 
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